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1 Introduction 



Abstract 

^ . In the framework of Heisenberg-Langevin theory the dynamical and sta- 

tistical effects arising from the linear interaction of two nondegenerate down- 
conversion processes are investigated. Using the strong-pumping approxima- 
tion the analytical solution of equations of motion is calculated. The phenom- 
' ena reminiscent of Zeno and anti-Zeno effects are examined. The possibility 

, of phase-controlled and mismatch-controlled switching is illustrated. 

o 

I Optical parametric processes yield a wide variety of optical phenomena. It 

■ is not surprising that many new phenomena will arise if a parametric process is 
^ I coupled to another one or to a different optical process. For instance, the super- 
O^' position of signal photons originating from two down-convertors with aligned idler 

beams leads to nontrivial quantum interference effects [Q. Parametric process cou- 
, , pled via Kerr interaction to an auxiliary mode, exhibiting quantum Zeno effect is 

rN I another nice example Many such composite systems (usually called nonlin- 

■ ear couplers) has thoroughly been studied in the literature. All-optical switching 
in the assymetric nonlinear coupler operating by the second-harmonic generation 
has been investigated in Q and its non-classical behaviour has been discussed in 
1^, ^. The quantum dynamics and statistics of the symmetric coupler containing 
two second-harmonic processes have been examined in |Q . The coupler composed of 
one linear waveguide and one nonlinear waveguide operating by the down-coversion 
process has been investigated in from the point of view of all optical switching. 
The occurrence of quantum Zeno and anti-Zeno effects in a similar device has been 
reported in ||l0|] . Amplitude behaviour of two linearly coupled down-conversion pro- 
cesses has been studied in [|j . Short-length analysis of this device has been given 
in §. ^ 

In this paper we deal with interesting phenomena arising as a consequence of 
linear interaction between beams propagating through the symmetric nonlinear cou- 
pler, which is composed of two nonlinear waveguides based on the down-conversion 
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Figure 1: Sketch of symmetric nonlinear coupler formed from two nonlinear wave- 
guides with susceptibility x''^^- The interacting beams are described by annihilation 
operators; L is the interaction length. 

processes. In fact this arrangement can be looked at as a continuous version of fa- 
mous Mandel's experiment , involving real physical interaction between the two 
down-conversion processes. In Section the equations of motion are derived and 
their analytical solutions are given. Sections || and |^ are devoted to the study of 
quantum dynamics and statistics of the coupler. Its non-classical properties are 
discussed in Section 



2 Equations of motion and their solution 

The coupler which is investigated in this article is composed of two nonlinear 
waveguides operating by the down-conversion processes in a directional arrangement 
(see Fig. The linear energy exchange by means of evanescent waves between 
pump, signal and idler beams is considered. The nonlinear media are assumed to 
be lossy. If such a system is far from resonance, the effective description involving 
only the field variables is adequate The effective momentum operator then 

reads 

G — Gi + G2 + Gres. + Gres.-syst. + Gint., (1) 



where 



h kjd^^aj + h (Viap.a^^ CLj +h.c^ for i=l,2, 

j=P,.,S„I, 

2 

Gres. ^ f^Yl H^y^l-^'J' 

i=l j=Pi,SiJi I 

2 



Gres.—syst. ^ 



Y 



i=l j=Pi,Si,Ii I 

Gint. = ?i(Kpapia^^ -I- Ksasi o^^ -I- K/a/ia|^ + h.c), (2) 

where aj(dj), j — Pi, P2, Si, S2, Ii,l2 are annihilation (creation) operators of pump, 
signal, and idler modes. Corresponding wavevectors along the z-axis of propagation 
are fcpj^, fcp^, ks^, ks2, ^/i and kj^. Linear coupling constants between pump, signal 
and idler modes are denoted Kp, ks and kj. Nonlinear coupling constants are 
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denoted as Fi and Each mode j is coupled via linear coupling constant kij to 
the Z-th reservoir mode characterized by annihilation (creation) operators bij{bj^) 
and wavevector kij along z-axis of propagation. The symbol h denotes the reduced 
Planck constant and h.c. represents Hermitian conjugate terms. 

The model represented by the momentum operator is symmetric both under 
the exchange 1^2, Kp <-» K*p, ks ^ '^gj <-> k*j and under the exchange 
S I. Since the dynamical behaviour of the system is completely determined by 
its momentum operator, the symmetries are conserved during evolution. This is 
convenient because it is not necessary to write down all calculated quantities, the 
rest being simply obtained by the above mentioned exchanges. 

Substituting (|l|) into the Heisenberg equations of motion (ih j^d = [G, a]), intro- 
ducing slowly varying operators Aj(z) — dj{z)exp{—ikjz) and applying the Wigner- 
Weisskopf approximation [ p^ , we arrive at the following Heisenberg-Langevin equa- 
tions of motion 



dAf 



dz 
dAs, 
dz 



-JPiAp^ + iK*pAp2exp{-iAkpz) + iriAsiAi^exp{-iAliz) + Lp^ (z), 
-75i^Si + iKgAs^expi-iAksz) + iTiAp^A\^exp{iAliz) + Lsjz), 



(3) 



where Afc^ = kk-^ — fc^^, k — P, S, I are linear mismatches, Ah — kp. — ks^ — kj., 
i = 1,2 are nonlinear mismatches, 7j, j = Pi, P2, 5*1, (S'2, /i, I2 are damping contants 
and the Langevin forces Lj{z) are assumed to be Markoffian 

= (i](z))-(L,(z)Lfe(z'))=0, 
(Lj(z)ifc(z')) = 2jj{ndj)Sjk6{z ~ z'), 

{L,{z)Ll{z')) = 2^,{{naj) + l)5,k5{z-z'). (4) 

Here angle brackets denote the averaging over the reservoirs, {udj) is one- mode 
mean photon number of the j'-th reservoir, 5jk is the Kronecker symbol and 5{z) is 
the Dirac delta function. It is useful to introduce the auxiliary quantities 



where 



and 



Ks^ = 7s, - lAKs^ , = lu + i^Ki^ , z = 1, 2, (5) 



Afc ± Aks Ak ± Akj 
= ^ : Aif/,^, = (6) 



2 

Afc 



1 ^ 

7^Y.^ks,+ki^-kp^). (7) 



2 

1=1 

The mismatch contains wavevectors of all modes and thus characterizes the 
overall phase mismatch. This important quantity will be called global mismatch in 
the following. 

If we assume the pump modes Pi, P2 are stimulated by the classical strong 
coherent fields 

ip,(z) (p,exp{-iAkpz/2), Ap^{z) -> ^P2exp(iAfcpz/2), (8) 

the system of equations of motion, represented by (^, splits into two independent 
sets. The first one corresponds to {^5^ , A5'2, , } operators and the second 
one corresponds to their adjoints. In what follows we will confine ourselves to the 
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first set. The special choice of the phases of classical amplitudes (H) leads, after the 
substitutions 

AsAz) = Cs,iz)eM~^^Ks,z), A\^{z) = C]^ (z)exp(^Ai^/,z), (9) 
to the system of linear differential equations with constant coefficients of the form 

= ^Ki,C\^~iKiC\^-iG\Cs,+t\{z), (10) 

where 

^sA^) = ^Si(2)exp(iAii'5iz), t\Az) = (z)exp(-iAis:/,z) 

are modified Langevin forces, and G\ = Ti^p^ , G2 — T2^P2 are rescaled nonlinear 
coupling constants. 

The system of Eqs. (^0|) can be solved using the Laplace transformation method 
and method of variation of constants. Returning to the operators Aj, the solution 
can be written in the following matrix form 

A(z) = M(z)[X(z)A(0) + R(z)], (11) 

where we have introduced the vector [( )^ means the transposition] 

A(z) = {AsAz),AsAz)JUz)Jliz)f, (12) 
the vector of reservoir contribution 

R(z) = iRs,{z),RsA^),RiAz),RiMf. (13) 
the diagonal matrix of mismatches 

M(z) = diag(exp(-iAii:siz) , exp(— iA_fi'52z), exp(iA_ft'/jz), exp(iA_ft'/2z)) (14) 

and the matrix of coefhcients 

4 

Xij{z) ^^{Ak)ijex.ip{Xkz) for i,j = l,..,4, (15) 

k=l 

where 

4 

Ak^ill {Xk-Xi)]-\Xla + \lh + \kC + d), k=l,..,A. (16) 

i=^k=l 

Four-dimensional matrices b, c, d can be found in Appendix |^ and a is unity 
matrix. The quantities A^, k ~ 1,..,4 in ( p^ ) and ( p^ ) are single roots of the 
polynomial 

A = x"^ + ax^ + bx"^ + cx + d, (17) 

with coefficients 

a = 7Si + 7S2 + 7/1 + 7/2 : 

b = Ls + Li + Li+L2+Ks,Ki,+Ki,Ks,, 

c = LsKi, + LiKs, + LiKs, + L2K1, 

Mks^'Ki, + \ki\''Ks, - \Gi\^Kj, - \G2\^Ks,, 
d = Ks,Kj,Ks,Ki, + \ks\^Ki,Kj, + \kj\^Ks,Ks,~\Gi\^Ks,Ki, 

-\G2\^Ks,Ki, + \kski -GIG2\^ (18) 
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where 



L,=Kj,K,, + \k,\\ j^S,I, L.^Ks^Ki^-lG^f, i = 1,2. 

If the roots of polynomial ([l7| ) are multiple, we can obtain the solution using the 
same methods. 



3 Quantum dynamics 

To investigate the dynamical behaviour of the coupler, we have to find roots of the 
characteristic polynomial ([TtI). 

If the damping is neglected and perfect phase matching is assumed, the poly- 
nomial is quadratic in and its roots are easy to find. In this case we can obtain 
periodical solution, exponentially amplifying solution or a combination of these two 
depending on the parameters of the process ]l3t . 

If either the losses are included and all phase mismatches are zero, or losses 
are neglected and phase mismatches are retained, we arrive at the fourth-order 
polynomial with real coefficients. The roots can be found using the Cardan formulae; 
unfortunately they are of complicated form and it is more convenient to solve for 
the roots numerically. 

In the most general case we need to solve the fourth-order equation with complex 
coefficients, a task, which can only be performed with the help of a computer. 
However, there are certain physically realizable regimes 

7Si = 7S2 = 7s, 7/i = 7/2 = 7/, 

Aks = Akj^O, Gi = G2^^, (19) 

ki\ks\ 

for which the general polynomial factorizes into two second-order ones. Their 
roots are 



-['ys + -yi + - |ks|)] ± Vbs -ii + iilKsl + \ki\ + Ak)]'^ + 4\Gi\ 
^1,2 = t: 



-[-ys 4- 77 + i{\fis\ - Ik/I)] ± ^y[Js - 7/ + i{\fis\ + - AA;)]^ 4- 410112 

^ 2 ■ 

(20) 

Neglecting the damping (75 = 7/ = 0), we can examine the influence of the global 
mismatch Afc to the dynamics of the coupler: 

1. If Afc e (-2|Gi|,2|Gi|), then 

(a) for \ks\ + G (0, |AA;| -I- 2|Gi|) all roots have the form = aj + ibj 
with nonzero real and imaginary parts aj and bj 

(b) for \ks\ + e (|AA:| + 2|Gi|,-|-oo) all roots are purely imaginary. 

2. If A/ce (-cx),-2|Gi|) U (2|Gi|,-|-oo), then 

(a) for |k5| + e (|AA:| — 2|Gi|, |Afc| + 2|Gi|) all roots are the same as in 
1(a) 

(b) for \ks\ + \ki\ G (0,|Afc| - 2|Gi|) U (|Afc| -f- 2|Gi|, -t-00) all roots are the 
same as in 1(b). 

Assuming the symmetrical linear coupling, \ks\ = the imaginary parts bj 

of Xj in cases 1(a) and 2(a) vanish and all roots acquire real values. In what 
follows, if Qj / {aj = 0), we will say that the coupler operates in the hyperbolic 
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\ks\ + |k/ 




Figure 2: Dependence of the character of dynamics of the coupler on global mis- 
match Afc and sum of linear coupling constants | K5 1 + 1 k/ | . Hatched area corresponds 
to the hyperbolic regime. The rest corresponds to the elliptic regime. 



(elliptic) regime. A more instructive demonstration of regions of different dynamical 
behaviour is given in Fig. |^. Notice first that Afc and its counterpart |ks|-|-|k/| affect 
the process in symmetrical ways. Now let us look closely at the interplay between 
the linear coupling and global mismatch. If |Afc| < 2|Gi| and the coupling strength 
is gradually increased, one can see that the coupler crosses the border between the 
hyperbolic regime (hatched area) and elliptic regime when a certain value of the 
linear coupling strength |k5| + is attained. More interestingly, for large phase 
mismatch |Afc| 3> SjCij (even such that almost no energy is converted from the 
pump mode to the signal and idler modes), the system moves up along the vertical 
line Afc = const, in Fig. ||, and it enters the region of instability characterized by 
the domination of the down-conversion part of the evolution when |k5| -I- Ik/I = 
Afc — 2|Gi|. For even stronger linear coupling the coupler leaves the region of 
instability again (crossing the line |ks| + |k/| = Afc-|-2|Gi|), the oscillatory character 
of the evolution is restored and pump photons gradually cease to decay. Interpreting 
(somewhat loosely) the linear coupling as a kind of continuous measuring process, we 
can look at the just described behaviour as being a manifestation of the well-known 
Zeno or anti-Zeno effects |j. Also here the strong influence of the "measuring 
apparatus" leads to the hindering of the decay of the originally unstable system. On 
the contrary under certain conditions (here nonzero phase mismatch Afc), the decay 
of the unstable system can be enhanced by a frequent (here continuous) monitoring 
of the unstable system 10 . Our Fig. || clearly shows the competition between 



these two opposite tendencies. 



4 Quantum statistics 

The quantum-statistical properties of the coupler are best studied employing the 
normal characteristic function containing complete statistical information about the 
system. The model represented by momentum operator (^ together with the lin- 
earization procedure (^) lead to the Gaussian characteristic function corresponding 
to the generalized superposition of coherent fields and quantum noise jl^ 



C^({/3,},z) = exp<^^ 



4 r , 
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4 

k=l,j<k 

+ (/3,e*(z)-c.c.)]}, (21) 

where the foUowing identification 5i = 1, 5*2 = 2, /i = 3, /2 = 4 has been done. The 
complex amphtudes ^i{z), ^2(2), Csiz), diz) are mean values of operators ^5^(2), 
(2:), Ai-^{z), Ai^{z), c.c. means the complex conjugated terms and 

B,{z) = (Ai](z)Ai,(z)), Q(z)^((Ai,(z))2), 
D,k{z) = (Ai,(z)Aife(z)}, = -(Aij(z)Aifc(z)) (22) 

for j.k — 5*1 , 5*2 , /i , /2 are noise functions. The complicated explicit expressions 
of the noise functions are given in Appendix The quantities Cj = Cj{0) and 
Bj = Bj (0) + 1 corresponding to the input beams are expressed under the condition 
of independence of incident beams in the form 

2 1 
Bj = cosh {r,j) + {uchj), Cj = -exp{i6j)smh{2rj), (23) 

where rj and 9j, j = 1,2,3,4 are squeeze parameters and phases of the incident 
beams and (rtchj) represents the mean number of external noise photons in the j-th 
mode. 

Assuming the unsqueezed input fields {vj = 0) the explicit expressions of noise 
functions (see Appendix ^) lead to the following identities 

CsAz) = CsAz) = CiAz) ^ Ci,{z) ^ 0, 
Ds.sAz) = Ds,i,{z) = Ds,iAz)^Q- (24) 

The quantum-statistical properties of single and compound modes can be quantified 
by means of many statistical quantities. From these we will use the principal squeeze 
variance X{z) [principal squeezing occurs if A < 1(2) for single (compound) mode], 

r - 1 ^ 

quadrature variances ( A|(z) ) [quadrature squeezing occurs if ((Ag)^) < 1(2) or 
((Ap)^) < 1(2) for single (compound) mode], normal reduced factorial moments of 
the integrated intensity jj^^^^^^^l — 1 [they are negative for non-classical states, neg- 
ative second moment reflects the sub-Poissonian photon statistics] and the photon 
number distribution p(n, z) [quantum oscillations in p(n, z) indicate the presence of 
state having no classical analogy]. 

Adopting the standard definitions of the above mentioned quantities for single 
mode [0 and using (p^), it is straightforward to show that single modes do not 
exhibit any interesting behaviour. In particular non-classical light cannot develop 
from coherent inputs in single modes. It can arise only as a result of quantum 
correlations of modes. 

In the case of the compound mode (j,j) the principal squeeze variance Aij(z), 
quadrature variances ([Agij(z)]^), ([Apij(z)]^) and variance of the integrated inten- 
sity ([AWy (z)]^) are defined as follows 

Kj{z) = 2{1 + B,{z) + B,{z) - 2Re[D,,{z)] - \a{z) + Q(z) -I- 2^,(^)1}, (25) 

((Ap2) =2{l + B,(z) + B,(z)-2Re[A,(2)]±Re[C,(z) + C,(z) + 2Aj(2)]} (26) 
and 

([AW^Az)]') = {[AW.{z)r) + {[AW,{z)r) + 2{AWdz)AW,iz)), (27) 
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where 



{[m^{^)f) = Bj{z) + + 2B,{z)Uz)? + 2Re[C,(z)C' W] (28) 



If the correlation function (|2^) is negative, we say that the corresponding modes 
are anti-correlated. 

The definitions of sum photon number distribution and k-ih moment (W^j(z)) 
are rather complex and can be found in |^|. 

5 Discussion of results 

As we have already mentioned above, analytical expressions of required quantities 
are only available under certain simplifying and restrictive assumptions. Even in 
those cases the expressions are of a complicated form and thus almost useless for 
qualitative discussions. Therefore we will employ numerical methods. The analyt- 
ical solutions, when available, may serve for checking the results of the numerical 
calculations. 

This section is devoted to the investigation of interesting phenomena arising from 
the linear coupling between two down-conversion processes. Each phenomenon is 
discussed in separate subsection. 

5.1 Quadrature switching 

It was reported in |^ that the symmetric coupler (|Gi| = IG2I) where only signal 
modes are linearly coupled (k/ = 0) behaves as follows. If mode 5*1 is squeezed in 
the given quadrature at the input, squeezing in a conjugated quadrature develops 
in mode 6*2. Taking into account also linear exchange between idler modes, we 
can observe a similar phenomenon in quadratures of compound mode {Si^Ii). Let 
us assume that both down-conversion processes are spontaneous, linear coupling 
constants are symmetric K5 = k/ and sufficiently strong. Changing now the phase 
ipp^ = arg^Pj of the pump mode P2 and leaving the phase ipp^ = arg^p^ of the 
pump mode Pi fixed, we can switch between quadratures at the output of mode 
{Si, Ii). Moreover, if the interaction length L is appropriatelly chosen, squeezing in 
the given quadrature can be transferred to the conjugated one in a continuous way 
(see Fig. I). 

5.2 Linear coupling can compensate wrong phases 

It is well known that for small interaction lengths z sub-Poissonian light can 
be generated in a nondegenerate down-conversion process in compound mode (5, /) 
provided that the process is stimulated (^SjC/ 0) and phases of incident beams 
fulfil the optimum phase condition arg(^5^/^p) = — ^. On the other hand, if 
either the process is spontaneous or the phase condition is strongly violated, this 
mode is super-Poissonian. Let us assume that the process in the first waveguide is 
stimulated by amplitudes S^Si i ^/i strongly violating the optimum phase condition 
(say arg(^5j^^/j^^p^) = ^) and the process in the second waveguide is spontaneous 
{^82 = C/2 — 0). Introducing the linear coupling between the waveguides, the 
modes (5*1, /i) and (82,12) can exhibit an interesting non-classical behaviour. The 
linear coupling restores the optimum phase condition and sub-Poissonian light is 



and 



(AW,iz)AW,{z)) 



2Re[D,,{z)^*iz)qiz) - D^.izMz)^ (z)] 
-f|A,(2)P + |A,(^)p. 



(29) 
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Figure 3: Demonstration of phase-controlled switching between the quadrature vari- 
ances ([Ag(L)]^) (o) and ([Ap(L)]^) (*) (the curve denoted by (A) corresponds to 
the principal squeeze variance X{L)) of mode (S'i,/i); L = 1.2, Ti = r2 = 1, 
Ks^Ki ^ 2, Afc = Aks = Aki = 0, ^p, = ICpJ = 1, = ^i, = ^s. = = 0, 
7j = 0.2, {udj) = 10-2. 



generated in mode {Si,Ii), surprisingly, for larger z (see Fig. ^ (a)). Further, sub- 
Poissonian light is also generated in mode (5*2,-^2) for small z (Figures ^ (b) and 
H), a phenomenon, which cannot be explained as easy as in the previous case. To 
deepen insight into this phenomenon we will resort to the analytical results. If losses 
are neglected, all mismatches are zero, k = Kg = k/ is real, G = Gi — G2, and 
K > \G\, we can calculate, using the results of Section |^ and (|2g|), the correlation of 
fluctuations 

{AWsJz)AWj,{z)) = |G|V(zK(z) + 2K2|G||^sJI6J"(^yW 

xsm{(pp - Lpsi - (fill), (30) 

where ifp, ipsn ^ii are the phases of input coherent amplitudes = ^p-^ = ^p^, 
Csi , 61 and 



u{z) 
v{z) 



COS 


V2 {^^ - 






3 sin 




-\G?) 


z 


4 









\ sm 



V2(^2_|G|2)2 



^2{k^-\G\^)z 



(31) 



Since both the expressions on the right hand side (R.H.S.) of Eqs. ( |3l| ) are real, 
only the second term on the R.H.S. of Eq. ( ^ can be negative depending on the 
sign of the product u{z)v{z), and on the argument of sine function. Restricting 
ourselves to small 2, we can expand u{z) and v{z) up to the z^ around the origin, 
and approximate Eq. (30) by the expression 



{AWsM^Wi,{z)) « \G\'z' + 2K'\G\\is^UiAz'^H^p - Vs^ - Vh] 



(32) 



Note first, that anti-correlation can only arise for sufRciently strong k and for suffi- 
ciently large z. It attains its maximum value if ^fip — ^Si^Vii — ^^siCp^Si^h) ~ 
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Figure 4: The integrated intensity {W{z)) ( — ) and its reduced factorial moments 

g^jMI _ 1 for k = 2{*),k = 3 (o), fc = 4 (A), A: = 5 (o) for mode Ji) (a) and 
(52,/2) (6): Ti = = 1, Kg = ki = 5, Afc = Afcs = Aki = 0, ^p, = ^p, = -i, 
^Si = 6i = 1, Cs, = ii, = 0, 7, = 0.2, (rid,) = 
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Figure 5: The sum photon number distribution p{n, z) for mode (5*2, ^2); the pa- 
rameters of the process are the same as in Fig. 0. 



(see Fig ^ (b)). It is also evident from the second term of the R.H.S. of Eq. ( ^2| ) 
that the Unear interaction enables us to affect the anti-correlation in mode (6*2, -^2) 
via amplitudes S^Si 1 ^ii ■ 

Repeating the arguments leading to the formula ( |30| ) for mode (Si, /i), we obtain 

{AWsAz)AWjAz)) = \G\'[2{\isA' + \^i.n + ^u'i^y{z) 

-2|G||^SillC/i|sin((/?P - V3si - ^h) 
xu(z)w(z) [u^{z) + |Gpw2(z)] . (33) 

Employing once more the expansion of (^Tj) around z — 0, we can see, that the 
second term on the R.H.S. of Eq. (33) is proportional to z and thus its sign is given 
only by the argument of sine function. Our choice of the initial phases then implies 
that the contribution of the second term is positive in this approximation. How- 
ever, the product u{z)v(z) alternates and its amplitude increases with increasing 
z, suppressing the quantum noise represented by the first term on the R.H.S. of 
Eq. (HI) and attaining the sub-Poissonian photon statistics for larger z. 



5.3 Cross mode 

Up to now we have separately discussed modes localized either in the first or in 
the second waveguide. This subsection will be devoted to the investigation of non- 
classical behaviour occuring in cross mode (6*1, /2). The investigation of the down- 
conversion with strong pumping led to the conclusion that signal mode S and idler 
mode / do not exhibit any non-classical behaviour irrespectively of the fact, if they 
are spontaneous or stimulated by the coherent light Q. Obviously, mode {81,12) 
compounded of modes 5*1 and I2 originating from two independent down-conversion 
processes cannot provide a non-classical light either. However, introducing the linear 
interaction between the processes, we can observe squeezing of vacuum fluctuations 
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Figure 6: (a) The quadrature variances {[Aq{z)]^) (o), {[Ap{z)]^) (*) and principal 
squeeze variance A(z) (A) for mode (81,12), (b) integrated intensity {W{z)) {- ) 
and its reduced factorial moments jj^^^^^^^l — 1 for = 2 (*), fc = 3 (o), & = 4 (A), 
fc = 5 (o) for mode (S'1,/2); Ti = r2 = 1, Kg = k/ = 4, AA; = Aks = Akj = 0, 
Cpi = = i, = = 1, ^s. = = 0, 7,- = 0.2, {udj) = 10-2. 
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and sub-Poissonian photon statistics simultaneously in this mode (see Figure 
To discover the origin of these phenomena we will employ once more the analytical 
solution. In the spirit of the derivation of Eq. ( |30| ) we can calculate the cross- 
correlation function 

Ds,iAz) = ->i\G\v''{z), (34) 

indicating, that linear exchange introduces the correlation between modes 5*1 and 
12- This correlation reduces both vacuum fluctuations (see (|25|)) and fluctuations 
of photon number. To make the latter more clear, we can derive the following 
cross-correlation function up to 

(AW^Si(^)AW^/,(z))«-2«;2|G||^sJIC/J^'sin(^p-^s, (35) 

It is worth noting that the effect of noise reduction can be enhanced by increasing 
the amplitudes and ^/^ . 



5.4 Mismatch- controlled switching 

Before discussing the last phenomenon we would like to mention several general 
remarks concerning the all-optical switching. This will enlight the motivation of 
the following discussion. Recent theoretical investigation of couplers has led to an 
interesting conclusion. Not only can they serve as a passive optical switchers, but 
they also provide the active control of the output of a particular waveguide by 
means of the input of the other one. There are at least two ways how to actively 
control the output beams. First, the coupling lenght of the coupler can be adjusted 
by changing the intensity of the strong classical input field Second, the 

phase-controlled distribution of the quantum noise in couplers can be realized p5[ 



(see also Subection 5.1). There is, however, one more possibility how to affect 
the properties of the outgoing beams. Inspection of Fig. || reveals that one can 
change the dynamical behaviour of the beams by means of the global mismatch Afc. 
Moreover, due to its global character (it contains all wavevectors) , we can control 
one mode by means of another one, even though they directly do not interact. The 
following arrangement can illustrate this. Let us assume that both processes are 
spontaneous, nonzero wavevectors fcs^, kg^, kj-^, kj^ and fcp^ are chosen in such a 
way, that they satisfy the matching conditions {Akg — Akj — Ali — 0) and linear 
interaction is in operation. Now the increase of the z-th component kp-^ of the 
wavevector of mode P2 entails the inhibition of the decay of the pump photons in the 
first waveguide (see Figure 0). This phenomenon is easy to explain based on Fig. ||. 
Initially, the parameters of the coupler are such that the down-convertion part of 
the evolution is dominant (we are inside the hatched area). The global mismatch 
Afc decreases with increasing fcp^ and the linear part of evolution grows dominant. 
Interpreting once again the linear interaction as a sort of continuous measurement, 
this effect can be looked at as a complementary effect to the Zeno or anti-Zeno- 
like effects described in Sec. ^. Unlike in Sec. || where initial condition (the value 
of Afc) was kept constant and the strength of "measurement" was changed, here 
the strength of the "measurement" is the fixed quantity and the initial condition 
is continuously varied. In this way the influence of the "measurement" results in 
the speeding up the down-conversion for small values kp^ and slowing down the 
down-conversion for larger fcp^ . This corresponds to a transition from the anti-Zeno 
to Zeno regime. 

It is also worth noting, that the integrated intensity of mode (82,12) depends 
on fcpj in the same way as mode (S'1,/1) does (see Fig. ^. This can be understood 
as follows. At the beginning (when fcp^ — 0) the process in the first waveguide is 
perfectly matched (A/i = 0) and the process in the second waveguide is strongly 
mismatched {AI2 ^ 0). The linear interaction, however, symmetrizes the device in 
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Figure 7: The dependence of the integrated intensity {W{z)) of ( — ) and 

(S'2, 12) (*) modes on the wavevector along the z-axis of propagation kp^ of pump 
mode P2; L = 1.5, Fi = r2 = 1, ks = ki = 2.5, ks^ = ks2 = 6, kj-^ = kj^ = 4, 



kp, = 10, Cp, = Cp, = 1, Cs, = = Cs, = 62 = 0, 7j = 0.2, {udj) = 10-2. 



the way that it partiaUy mismatches the first process and partially compensates the 
mismatch in the second waveguide at the same time. 

6 Conclusion 

The quantum dynamics and statistics of the symmetric nonlinear coupler operating 
by down-conversion process have been investigated. In a framework of strong pump- 
ing approximation we have solved analytically the Heisenberg-Langevin equations. 
The manifestation of Zeno and anti-Zeno effects has been demonstrated based on the 
analytical solution. The non-classical behaviour of beams involved has been studied 
based on numerical calculations. The phase-controlled redistribution of quantum 
noise between the quadratures can be achieved in mode {Si , /i ) . The possibility 
of generation of sub-Poissonian hght in modes {Si,Ii) and {82,12) caused by the 
linear interaction of two super-Poissonian lights has been shown. Light exhibiting 
simultaneous squeezing of vacuum fluctuations and sub-Poissonian photon statis- 
tics can be obtained in cross mode (5*1, 12). The inhibition of the decay process in 
the first waveguide owing to the nonlinear matching of the second process has been 
observed. All these phenomena were shown to be robust against the presence of 
weak damping. 



A Matrices b, c, d of Eq. 



b = 




iK*s 

Ksi + Ki, + 



id 






iG2 

— IKI 











—ik'i 



Ks, + + Ki, 
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d = 



Ki^ +L1 + L2 iK*s {Ki^ + Ki^ )_ iGi {Ks2 + Ki^ ) 

-iGliKs:, + Ki^) - KiG*2 Ks^Ki^ + Ls + U 

KsGl - K'iGl -iGl{Ks^ + Ki,) ^in*r{Ks, + Ks,) 

iG2{Ks,+Ki,) 
-iKi{Ksi + KS2) 
Ki^Ks2 +Ls + Li 



L2Ki,+\ki\^Ks, 

insLi — in*iGiG2 



iK*gLi — iKiGiG2 iGiL2 + iK*gK*jG2 

LiKl, + [kiI'^Ks, k*jG2Ks, - ksGiKl, 

-iGlL2 - iKsKiGl KsGlKj^ - K1G2KS1 L2Ksi + \KsfKi2 

Gl —iK*iLs + insGiG'2 

-K*sG2Kl^ + K1G1KS2 

iG2Li + iKsniGi 

—IKlLs + iK*sGlG2 

LiKs2 + iKsl'^-ftT/i 



B Noise functions 



Bs,{z) = {AAl^{z)AAs,{z))=^{\X,j\^B, + 2^j{ndj)xij) 



+ E(27j+2Xy+2-|Xi,f), 



Cs, {z) = ((Ais. {z)f) = E i^ii^j + exp(-2»Ai^5.^), 
DsiS2{z) 



{AAsAz)AAs2iz)) = ^ + Xy+2X2,+2q+2) 

xexp(-iAfc), 

4 



Ds,i,{z) = {AAs,{z)AAi,{z)) = 



exp[-i{AKs,+AKi,)z], 



Ds,i2{z) = {AAs,{z)^Ai2{z)) 



+E(27.-xii-^y+2n-+2) 



eM-i{^Ks,+AKi2)z], 



Ds,s2(.z) = -{AAliz)AAs2{z)) ^ 



J2{XljX,^B^ + 2^^{n,^)xfj) 



+ ^{2^,+,xf+2j+2-XljX,,) 



exp(— iAfcs^), 
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2 

Ds,I, {z) = -(A4^ (z) Ai,, (z)) - - ^ (^r,^3,q + XiV2^3V2Q+2) 

J = l 

xexp[i(Ais:s, - AX/Jz], 

2 

Ds,I, [Z) = -(A4^ (z) Ai,, (Z)) = - ^ (^r,^4*,q + ^r,+2^1,+2Q+2) 

3 = 1 

xexp[i(Aifs, - AX.Jz], (36) 



where = Xij (z) are defined in Eq. ( p^ ) and 

Xij = Xij{z)^ j \Xij{z - z')\'^dz', 
Jo 

Xfi = xfiiz) ^ f X,,{z ~ z')XUz - z')dz' . 
Jo 

The rest of the noise functions can be obtained using the symmetry of the model. 
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